Magnetism and Mott Transition: A Slave-rotor Study 
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Motivated by the debate of spin-density-wave (SDW) versus local-moment (LM) picture in the 
iron-based superconducting (FeSC) materials, we consider a two-band orbital-symmetric Hubbard 
model in which there is robust Fermi surface nesting at (7r,0). We obtain the phase diagram 
of such system by a mean-field slave-rotor approach, in which the Fermi surface nesting and the 
SDW order are explicitly taken into account via a natural separation of scale between the Hund's 
coupling and the Coulomb interaction. We find that for a sizable range of Hund's coupling the Mott 
transition acquires a strong first-order character, but there also exists a small range of stronger 
Hund's coupling in which an enhancement of magnetization can be observed on the SDW side. We 
interpret the former scenario as one in which a sharp distinction can be drawn between LM and the 
SDW picture, and the latter scenario as one in which signs of LM physics begin to develop in the 
metallic phase. It is tempting to suggest that some FeSC materials are in the vicinity of the latter 
scenario. 
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I. INTRODUCTION 

Mott transitions in multiband scenarios have recently 
received increased attention in the condensed matter 
community, partly owing to the interests in materials re- 
lated to the iron-based superconducting (FeSC) materi- 
als, whose parent compounds are believed to be in close 
proximity to Mott transitions^^ Indeed, while most par- 
ent compounds of the FeSCs are poor metals, insulating 
behaviors have been observed in^ (Tl,K)Fea;Se2 and in^ 
La202Fe20(Se,S)2. Given the sizable magnetic moment 
of approximately 0.8-1.0 /is in the parent compounds of 
the "122" family FeSCs^'"— a persistent debate in the 
field has been whether the magnetism is best described 
by an itinerant spin-density- wave (SDW) nestin g^^i^^ or a 
local-moment (LM)^ii^ picture. Theoretically, ah initio 
LSDA calculationi^ is able to obtain the experimentally 
observed magnetic ordering pattern in these compounds, 
but the magnetic moment is overestimated to be approx- 
imately 2.0 suggesting that a weak-coupling SDW 
picture alone may not be adequate in describing these 
materials. 

With the above considerations, it is beneficial to con- 
sider how magnetism and Fermi surface nesting affect 
the Mott transition, and vice versa, in a multiband sce- 
nario. Of particular interest is the question of whether 
there are any signs of local- moment physics on the metal- 
lic side of the phase diagram. However, while multi- 
band Mott transitions have previously been studied via 
dynamical mean-field theorjii^^— and slave-spin mean- 
fieldfi^ many of these studies have been focused on para- 
magnetic states. 

In contrast, in this paper we present a slave-rotor— 
study of an orbital-symmetric multiband Mott transition 
in which, utilizing a natural separation of scale between 
the orbital-symmetric Coulomb repulsion and the Hund's 
coupling, the SDW Fermi surface nesting has been explic- 
itly taken into account. We find that for a sizable range of 



Hund's coupling, the existence of (nearly) nesting Fermi 
surfaces causes the Mott transition to acquire a strong 
first-order character, in which the (staggered) magneti- 
zation jumps across the phase boundary. However, as 
the strength of Hund's coupling further increases, the 
first-order transition becomes weaker and an enhance- 
ment of magnetization can be observed on the metallic 
side. We interpret the former scenario as one in which 
a sharp distinction can be drawn between the itinerant 
SDW nesting and the LM picture, and the latter scenario 
as one in which signs of LM physics begin to develop in 
the metallic phase. 

It is worth noting that in the region of parameter 
space in which the staggered magnetization is enhanced, 
the renormalized hopping parameters also develop a dis- 
cernible anisotropy between the x and y nearest-neighbor 
bonds, in which the antiferromagnetic x direction shows a 
larger renormalized hopping than the ferromagnetic y di- 
rection. Such anisotropy is insignificant in other regions 
of the parameter space we considered. In the context of 
FeSC, such anisotropy may be associated and is in agree- 
ment with the observed resistivity anisotropy in detwined 
Ba(Fei_2;CO:c)2As2 samples;^ in which the in-plane re- 
sistivity in the antiferromagnetic direction is found to be 
smaller than that in the ferromagnetic direction — a result 
difficult to reconcile with a simple Landau Fermi liquid 
picture)^ 

We also remark that the renormalized bandwidth we 
obtained in this region of enhanced magnetization is 
about 2.5-5 times narrower than that of the correspond- 
ing noninteracting tight-binding model — values that are 
close to, albeit larger, than those reported for some FeSC 
materials i^^iSi Similarly, the staggered magnetization in 
such region is around 1.3-1.5 /is, close to but slightly 
larger than the values observed in the FeSC materials in 
the "122" family. Together with the observed resistiv- 
ity anisotropy, it is tempting to suggest that some FeSC 
materials in the "122" family are close to this region of 
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enhanced magnetization that we associated with remnant 
of LM physics in the metaUic phase. 



II. THE MULTIBAND HUBBARD MODEL AND 
THE SLAVE-ROTOR FORMULATION 

A muhiband Mott transition can be described by a 
muhiband Hubbard model, which can be written as the 
sum of a noninteracting hopping Hamiltonian and an on- 
site interaction term: 

H = Hhop + Hint , (1) 

Hhop = Yit'i^claa'^jba + h.C.) ^ M XI ^laa^iaa , 

i,j a,b i a 

(2) 

Hint = {y'Y^AaaAaa''^ia,iy'Ciaa 



I a 

U -2Ji 



Et t 



a^b 



(3) 



where i,j label lattice sites, a,b label orbitals, and a, a' 
label spins. Note that the sum over repeated spin indices 
is assumed in the above equations and will be assumed 
throughout this paper. Note also that we assume U to 
be orbital independent, as would be case the when all 
orbitals are atomic. In that limit, it can also be shown^^ 
that Ji = J2 ^ J3, which we shall assume henceforth and 
thus drop the subscripts. Henceforth we shall also refer to 
the U term as the Coulomb interaction and the J term as 
the Hund's coupling, even though the latter includes the 
contributions from co-hopping and inter-orbital Coulomb 
interactions. 

For concreteness we take i/hop to be the orbital- 
symmetric two-band Hamiltonian introduced by Ran et 
al. in Ref. [26l . whose hopping parameters and Fermi sur- 
face topology are shown in Fig. [TJ Such Hamiltonian 
is believed to capture certain essential physics in FeSCs, 
and in such context the two orbitals can be understood as 
the dxz and dyz orbital of the Fe site, with the X, Y axes 
oriented in the Fe layer and at 45° with the Fe-Fe bond 
direction. As is evident from Fig. l(b)[ at half-filling 
(i.e., two electrons per site) there is strong Fermi surface 
nesting in -ffhop, with nesting vector (tt, 0) and (0, tt). con- 
sistent with the SDW ordering in theSL^^ "1111" and^i^ 
"122" families FeSCs. However, as explained in Ref. [2^, 
for topological reasons the Fermi surfaces are not fully 
gapped even when interactions are included; hence there 
is a robust metallic SDW phase in such model. 

In the slave-rotor formalism^ the Hilbert space is en- 
larged by decomposing the electron operator Ciaa into 
an 0(2) rotor 6i and a fcrmionic spinon fiaa, such 





(b) 



FIG. 1: (a) The hopping parameters of the two-band hopping 
Hamiltonian Hhop in Ref. [2^; (b) The resulting Fermi surface 
for ti = l,t'i = 0.2, t2 = 1.7, and 4 = 0.3 at half-filling (i.e., 
two electrons per site), where blue represents hole pockets 
while red represents electron pockets. The choice of param- 
eters and filling fraction in (b) are assumed throughout this 
paper. 
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Note that in such decomposi- 



tion there is only one rotor per lattice site and the spin 
and orbital indices are carried solely by /. Such econ- 
omy in the slave-rotor formalism is made possible by 
the orbital-symmetric structure of the Hubbard Hamil- 
tonian Eq. [TJ particularly in the existence of a large 
orbital-symmetric Coulomb interaction U . In this en- 
larged Hilbert space, the physical subspace is given by 

the constraint Y.a Y.a [fLaf^aa- - 1/2) = Li, where Li 
is the angular momentum operator conjugate to 9i. 

Plugging the slave-rotor decomposition into the Hub- 
bard Hamiltonian Eq. [TJ and applying a mean-field ap- 
proximation to decouple the spinon / and the rotor 6, we 
obtain: 

H « Hmf = Hj: + He , (4) 

= E E(*?^/L/.b. + h.c.) - E 5^(M+/iO/L/.a. 



i,j a,b 

J 
2 



■ ( '^fiaafiba' fiba' fiacr 

i a^b 

+ fiaa fiba' fiair' fiba + fjaa fiaa' fiba'fibaj, (5) 



He = YP^^^""'-''^ + h-c-) + E ( + h^L^ ) ' 



(6) 



where hi are mean-field parameters introduced to en- 
force the constraint on average, f^j* = ifj'(e'*^^'~^^)), and 

= Y.a,b*i] (fiacr fibty)^ all of which are to be deter- 
mined self-consistently. In this paper we consider only 
translationally invariant solutions in which the rotor unit 
cell is unenlarged while the spinon unit cell is enlarged to 
include two lattice sites per cell. Moreover, since the un- 
derlying Hamiltonian is time- reversal invariant, we con- 
sider only solutions for which both J'ij and tf^ are real. 
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Furthermore, since the system is at half-fihing, the con- 
straint is satisfied by setting hi = 0. In this paper, the 
self-consistency for t^j and Jij arc solved by repeated 
iterations. 

Following Florens and Georges;^ a second mean-field 
approximation is applied to the rotor Hamiltonian Eq. |6j 
in which the rotor variable e**' is replaced by a complex 
bosonic field Xi subjected to the constraint \Xi\'^ = 1, 
enforced on average by a second mean-field parameter A. 
With such approximation the rotor Hamiltonian Eq.[6]rc- 
duces to a system of coupled simple harmonic oscillators. 
The electron quasiparticle weight in this formulation is 

given by Z = 1 - / ^ y^ 4(£x(k)+A) ' ^^ere £x (k) is the 
dispersion of the X boson in k space. The system is in 
the metallic phase when Z ^ 0, which happens when the 
X boson condenses at its energy minimum. 

For the spinon sector, we apply mean-field factoriza- 
tion to the four-fermion terms that appear in the spinon 
Hamiltonian Eq. [5] Specifically, given the Fermi sur- 
face nesting, we consider the following SDW mean-field 
Hamiltonian: 



i?sDw = Y.T.(^'fLf,b. + h.c.) - /i E E /L/.a. 

i,j a,b i a 

+ E(-i)^^ E ^'^'^b (/-t/^H - fliUi) , (7) 

i a,b 
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FIG. 2: The phase diagram of the half-filled orbital-symmetric 
two-band Hubbard Hamiltonian Eq. [U where J and U are 
given in units of ti (for comparison, note that the full band- 
width of the noninteracting system is VK = 12.8 ti). Here P 
stands for paramagnetic metal, SDW stands for spin density 
wave metal, LM stands for antiferromagnetic insulator, and 
SL stands for spin liquid. The insulating region is indicated 
by the gray shade, and first-order phase boundaries are indi- 
cated by the yellow shade. For an explanation of the "NNN 
decoupling" and "NN decoupling" lines, see Fig. [4] and the 
main text. 



III. RESULTS FROM SLAVE-ROTOR MEAN 
FIELD 



where denotes the x coordinate of site i, and Mat 
is a Hermitian matrix of mean-field parameters, to be 
determined for given t^j* by minimizing (Hf) with respect 
to the mean-field state obtained from i?sDW- Note that 
since Mab carries orbital indices, the possibility of an 
orbital symmetry broken SDW state is included in our 
mean-field ansatz. 

Since the -ffsDW breaks the symmetry between the x 
and y direction in the underlying Hamiltonian Eq. [l] the 
self-consistent mean-field parameters Sij and Uj can take 
different values along the x and y direction. We have 
also allowed for the possibility in which jTij and Uj take 
different values along the {x + y) and {x — y) directions, 
but found no manifestation of such cases in the parameter 
space we considered. 

It is well-noted that different parts of the electron two- 
body interactions arc treated differently in our formu- 
lation, with the Coulomb interaction U handled via a 
"strong-coupling" approach while the Hund's coupling J 
handled via a "weak-coupling" approach. While this may 
seem unsatisfactory at first sight, it may nonetheless cap- 
ture the qualitative aspects of the system when U is the 
dominant scale, which seems to hold true for 3d transi- 
tion metals^ and has been assumed in other theoretical 
calculations involving FeSCsi^i^i^i^. 



In this section we present the results of the mean- 
field calculations outlined in Sec. [iTj for the half-filled 
orbital-symmetric two-band Hubbard Hamiltonian Eq. [TJ 
without worrying about the validity of the approxima- 
tions made. We shall return to the issue of validity 
in Sec. IIVI To simplify notations, we measure both J 
and U in units of ti (for comparison, note that the full 
bandwidth of the noninteracting Hamiltonian Eq. [2] is 
W ^ 12.8 <i) and denote 0{,- = {X^X*) = (e^C'--"^)), 

such that = Qljtl^- We shall refer to Q/^- as the bond 
renormalization factor, t^j* as the renormalized hopping 
parameter, and the noninteracting part of Hf in Eq. [5] 
as the renormalized hopping Hamiltonian. We also de- 
note the orbital-dependent staggered magnetization as 

niab, i.e., mab = Ei(-l)*"(/ItM - fiaJibillJ^ (^ere 
N is the number of lattice sites), and decompose it as 
mab = \{mo5ab + "T-iTab + "^2T^b + m^.Tl^) , with being 
the Pauli sigma matrices. 

Our major result is the phase diagram shown in Fig. [21 
in which we vary J and U while keeping the bare hop- 
ping parameters t°ij fixed. To provide further insight 
into the phase diagram, we also plot in Fig. [3] the quasi- 
particle weight Z ^ the bond renormalization factor Q^- 
along the x, y, and diagonal {x + y) direction, and the 
orbital-diagonal staggered magnetization mg as function 
of U for four characteristic values of J . It should be 
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FIG. 3: The quasiparticle weight Z, the bond renormahzation 
factor Q{j along the x, y, and diagonal {x + y) direction, and 
the orbital-diagonal staggered magnetization mo as function 
of U for (a) J = 0.4, (b) J = 0.8, (c) J = 1.2 and (d) J = 2.2. 
The inset in (d) shows a blowup of the bond renormahzation 
factors in the region where magnetization is enhanced. 



noted that within the parameter range we considered, 
the orbital-dependent staggered magnetization niab is al- 
ways dominated by the orbital-diagonal component mn, 
in agreement with the weak-coupling results in Ref. l26l . 
In other words, no tendency towards an orbital-selective 
SDW state is found in our results. Because of this, hence- 
forth we shall consider exclusively mg, and refer to it 
simply as the staggered magnetization. 

As typified by Fig. |3(a)[ for J < 0.5 we find that the 
system undergoes a second-order Mott transition before 



KIXIXIX 
<IXIXI>< 

XXIXX 

(a) 




(c) 



FIG. 4: The real space picture of the renormalized hopping 
Hamiltonian on (a) the left of the "NN decoupling" line, 
(b) between the "NN decoupling" line and the "NNN decou- 
pling" line, and (c) the right of the "NNN decoupling" line in 
Fig. [2] Here the green dots denote the location of Fe atom, 
the violet (blue) lines denote non-zero renormalized nearest- 
neighbor hopping along the x (y) direction, the brown lines 
denote non-zero renormalized next-nearest-neighbor hopping, 
and the gray dotted lines denote zero renormalized hopping. 



developing any magnetism. The character of the Mott 
transition is essentially the same as the one considered 
by Florcns and Georges in Ref. [l^- As a result, for small 
values of J part of the parameter space supports a non- 
magnetic insulating state, i.e., a spin liquid. Interest- 
ingly, the onset of magnetism for this parameter range 
seems to coincide with a phase boundary (labeled as 
"NN decoupling" in Fig. [5]) in which the nearest- neighbor 
bonds renormalize to zero. 

As J increases to0.7<J<0.9, the onset of mag- 
netism begins to shift to the left of the Mott transi- 
tion. Moreover, the first-order boundary in which the 
nearest-neighbor bonds renormalize to zero starts to co- 
incide with the Mott transition. The combined transition 
retains a first-order character, and the staggered mag- 
netization exhibits a discontinuity across the transition. 
Furthermore, another phase boundary in which the next- 
nearest-neighbor bonds also renormalize to zero (labeled 
as "NNN decoupling" in Fig. [2]) can now be seen, whose 
critical value ?7^^^ decreases with increasing J. Note 
that the phase boundary of this next-nearest-neighbor 
decoupling traces well with the onset of a fully polarized 
insulating state. This situation is typified in Fig. |3(b)[ 
For clarity, we also illustrate the real-space picture of 
the nearest-neighbor and next-nearest-neighbor decou- 
pled state in Fig. |4l 

As J further increases to 1.0 < J < 2.0, the Mott 
transition becomes even stronger in first-order character 
and starts to coincide with the transition in which the 
next-ncarcst-neighbor bonds decouple. The phase tran- 
sition is now accompanied by a jump to full polarization 
across the phase boundary. This situation is typified in 
However, as J further increases to J > 2.1, 



Fig. 3(c) 



the first-order character of the Mott transition starts to 
soften up, and, as typified by Fig. |3(d)[ a small shoul- 
der in which the staggered magnetization is enhanced 
can now be observed immediately to the left of the Mott 
transition. 

Moreover, as can be seen in the inset of Fig. |3(d)[ an 
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(a) (b) (c) 

FIG. 5: Three patterns of renormalized bonds in a square 
lattice Hubbard model with non-mixing bands having only 
nearest-neighbor hopping. Here the green dots denote lattice 
sites, the blue line denote nonzero renormalized hopping, and 
the gray dotted lines denote zero renormalized hopping. It 
can be checked that within the slave-rotor mean field and on 
the insulating side (c) is energetically favorable to (b), which 
in turns is energetically favorable to (a). 



anisotropy of bond rcnornialization factors between the x 
and y nearest-neighbor bonds has also become discernible 
in this region where staggered magnetization is enhanced. 
In contrast, such anisotropy remains insignificant in all 
other parts of the parameter space we considered, includ- 
ing regions in the metallic phase in which the staggered 
magnetization is not enhanced. Note also that the bond 
renormalization factor is larger in the antiferromagnetic 
X direction than in the ferromagnetic y direction. 

IV. DISCUSSIONS AND CONCLUSIONS 

Since the results in Sec. IIIII rely on various mean-field 
approximations, one should be cautious about the results 
thus obtained. In particular: 

1. The critical value of U for the Mott transi- 
tion, obtained in the slave- rotor approach 
is smaller than that obtained in the slave-spin 
approachii^ This is a known discrepancy between 
the two methods and should vanish in the large- 

(number of orbitals) limit. For the case where 
N is finite the slave-spin estimate is believed to be 
more accurate<^ This, however, should not affect 
the qualitative statements made in this paper. 

2. In the slave-rotor formalism, all bonds that con- 
nect between the same sites are renormalized by 
the same factor, regardless of orbital character. 
This may lead to worries that some orbital selec- 
tive aspectsi^iii of the system are neglected. Such 
worries are partially alleviated by the stud}*i^ of 
a similar two-band system in the slave-spin for- 
malism, which allows for the possibility of spin- 
and-orbital dependent bond renormalizations and 
finds no evidence for an orbital selective state for 
J > 0. Furthermore, we reiterate that the possibil- 
ity of an orbital symmetry broken SDW state is in- 
cluded in our formulation via an orbital-dependent 
mean- field parameter Mab, but that we always find 



the staggered magnetization of the resulting state 
to be dominated by the orbital-diagonal compo- 
nent. Therefore, while orbital-selective scenarios 
have been proposed in the context of FeSCsj^i^i^i 
we have found no supporting evidence for such state 
in the two-band model we considered. 

3. As explained in Sec. |TT1 the distinct treatment for 
the Coulomb interaction and Hund's coupling is 
justified by the assumption that U is the domi- 
nant scale. Generally, spin fluctuation is believed to 
play a more important role when the inter-orbital 
Coulomb interaction U' = [/ — 2J is comparable 
to J, i.e., when U < 3 J. Since 

[/Mott decreases 

rapidly as J increases, the phase transition starts 
to enter the region where U < 3J around J « 2.5. 
Therefore, the upper part of the phase diagram in 
Fig.[5]may be stretching the limit of validity for the 
present approach. 

4. By treating J via a weak-coupling approach, we 
have neglected the strong-coupling aspect of the 
Hund's coupling. In particular, one should not ex- 
pect the to remain unchanged for small J 
as in Fig. [5J Instead, one should expect U^°** to 
drop immediately as J increase, as is foundi^ in the 
slave-spin approach. 

5. In similar spirit, the present approximation has ne- 
glected the weak-coupling aspects of U. In particu- 
lar, at the perturbative level both U and J are ex- 
pected to drive the Fermi surface nesting at (tt, 0). 
Hence, the paramagnetic metal and the spin-liquid 
region in Fig. [2] should, if exist at all, be greatly 
reduced in size. 

6. One peculiar feature of the mean-field result pre- 
sented in Sec. IIIII is the existence of phases in 
which certain bonds in the hopping Hamiltonian is 
renormalized to zero in the insulating phase. This 
seems to be a general feature of the slave-rotor 
method, since similar results is also found in a Hub- 
bard model with a different lattice geometry (and 
which the Hund's coupling is dropped) Indeed, 
the problem can already be seen in the simplest case 
of a square-lattice Hubbard model with non-mixing 
bands having only nearest-neighbor hoppings. In 
such case it can be checked that on the insulating 
side, within the slave-rotor mean field, a symme- 
try broken phase in which all bonds along a par- 
ticular direction is renormalized to zero is energet- 
ically favorable to a state in which all bonds are 
renormalized by the same amount. This rotation- 
symmetry broken state is, in turn, less energetically 
favorable to a translation-and-rotation-symmetry 
broken state in which the bonds renormalize to a 
dimer pattern (see Fig.[5]for illustration). Since the 
dimer pattern is known to be favorable in the large- 
N limit )^^'^^ the decoupled phase in Fig. [5] may be 
an artifact of the large- approximation. 
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In spite of the issues mentioned above, the present 
study may still shed light on the question we posed in 
the introduction of this paper: namely, whether there 
are any signs of local-moment physics on the metallic 
side of the phase diagram. For this, observe that C/^°" 
remains in the region where C/ > 3J for J in the range 
0-7 ^ J ^ 2.5, and as we approach the phase bound- 
ary from U < the bond renormalization remains 
modest for much of the parameter space, such that the 
renormalized bandwidth W is at least comparable to J. 
In such region the approximation we took in this paper 
may be valid. And for 0.7 ^ J ^ 2.1, as we approach 
the phase boundary from U < [/^°** there is no sign 
of magnetization enhancement (which, in this formula- 
tion, will be accompanied by a sharper-than-usual drop 
in Qj.;) until the first-order phase boundary is hit. Thus, 



our result is suggestive that in this range of J there is 
no sign of local-moment physics on the metallic side. In 
contrast, for 2.1 < J < 2.5 there is a sharper drop in 
Qlj as we approach the phase boundary, which suggests 
the possibility that in this range of J some remnant of 
local-moment physics can be found. 
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